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Application of Fokker-Planck equation to heavy quark transport in the evolving medium created 
in heavy ion collision is critically scrutinised. We realise that the approach introduces a moderate 
uncertainty in drag and diffusion coefficients culminating in huge ambiguity in the theoretical pre¬ 
diction of nuclear modification factor Raa . Quantitative estimation of the error is presented by 
considering recent developments in this field. 


The last thirty years have witnessed fervent activ¬ 
ity in the physics community aiming to understand 
de-confined quark matter created through high energy 
nucleus-nucleus collisions. Ideas started to percolate 
when studies of finite temperature quantum chromody¬ 
namics (QCD) [l":3) suggested the existence of a phase 
where the quark and gluonic degrees of freedom appear 
instead of being confined within hadrons. Another ad¬ 
vertisement that has been inviting further attention to 
heavy ion collisions (HIC) is the scope to study our early 
Universe from the convenience of a laboratory. But, one 
should be cautious because in early Universe the evolu¬ 
tion of quark gluon plasma (QGP) was quite slow com¬ 
pared to that in HIC ifzj. The time scale of evolution of 
QGP in early Universe is expected to be of the order of 
micro-seconds [s] whereas the QGP in a laboratory does 
not even persist for more than 10 — 12fm j|. It is, there¬ 
fore, a matter of debate how far one can infer about the 
early Universe from a study of quark gluon plasma in a 
laboratory where it is eventually subjected to a violent 
evolution. But this ambitious attempt may wait! First 
and foremost, we need to determine how well such a fast 
evolving medium can be understood with the tools at our 
disposal. 

Clear and convincing evidence for the formation of a 
strongly interacting medium was announced by experi¬ 
ments in 2010 Q. Theoretical and experimental projects 
have been on run to figure out various properties of 
this medium. To explore such a short-lived medium, 
one has to rely on certain internal probes . Ones like 
heavy quarks are born much before the creation of QGP 
and witness the whole of its evolution. Essentially, the 
heavy quarks lose energy inside the medium and this 
phenomenon of energy loss, the jet quenching, is a well 
established signal for the formation of such a medium 
[|j. Obviously, this energy loss relies on the transport 
of heavy quarks through the medium. With the help 
of transport coefficients, the change in initial momen¬ 
tum distribution, /)„, of heavy quarks can be computed 
by using Boltzmann transport equation or Fokker-Planck 
equation (FPE). Hence, the relative change in distribu¬ 
tion -p-, ff designating the final distribution of heavy 
quarks, yields the nuclear suppression factor, Raa, an 


experimentally measurable quantity. 

The FPE can be represented in the following form: 


dl = d_ 

dt dpt 




(1) 


where p characterises momentum of the heavy quark with 
a momentum distribution / . A and B are the drag and 
diffusion coefficients respectively at a fixed temperature 
of the medium. Although, in heavy ion collision it is still 
not beyond suspicion whether the medium gets enough 
time to thermalise, yet to make life easier one usually 
advocates for a local thermal equilibrium. The applica¬ 
tion of FPE to a static QGP is very straight forward for 
momentum independent transport coefficients (drag and 
diffusion) @ and the solution to eqjTJhas been deduced a 
number of times in this context. Some studies have been 
done by considering momentum dependent drag and dif¬ 
fusion as well 0 - 

The scenario in heavy ion collision is drastically differ¬ 
ent from that of a fixed temperature thermal bath be¬ 
cause the temperature of QGP falls off rapidly as the 
fireball expands. A clever way of incorporating this into 
the Fokker Planck equation has been adopted in recent 
literature for both momentum independent 

and momentum dependent transport coefficients. There, 
the drag A and diffusion B in eqJT] are replaced by time 
dependent functions owing to their implicit dependence 
on time through the temperature. This scheme tacitly 
introduces a time scale At within which the variation in 
temperature is so small that Fokker Planck equation for 
a fixed temperature bath can be safely applied. In princi¬ 
ple, At must tend to zero to justify the replacement of A, 
B by A(t), B(t). But, it is not permissible as the FPE is 
only valid at a time scale r which should be much bigger 
than the collision time [] and much smaller than the 
relaxation time of heavy quarks in medium. Hence 
At should at least be of the order of r. Furthermore, 
to guarantee collisions, At should be considered not less 


1 t ( H ® is the time for which the heavy quark comes into the ’ 

sphere of influence’ of the bath particle 
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than the heavy quark mean free time Tj v . So, we need to 
check whether these considerations are compatible with 
the present scenario of rapidly evolving QGP in heavy 
ion collisions. 

In centre of momentum (COM) frame, can not 
be less than the inverse of average transverse momen¬ 
tum transferred per collision. We go with the usual con¬ 
vention by considering a gluon dominated plasma and 
hence, the averaging may be performed with the help 
of the heavy quark (Q)-gluon(g) differential scattering 
cross-section, 


da 

dt 


27ror 


(i-a 2 m ) 2 


( 2 ) 


where a is the strong coupling. For a heavy quark with 
mass M, COM collision energy \/I, A 2 M = M 2 /s. s and 
t are the Mandelstam variables. 

Using eq. [2] the square of the average transverse mo¬ 
mentum transferred is derived to be, 


<«!) = 


sm 


D 


-In- 


s — 4 m 2 D 4 m 2 D 


( 3 ) 


Here mu = 


gT ^I (Ca + ^ with g = y/Anxa is the 

Debye mass at temperature T. Ca = 3 is the Casimir 
factor and Nf = 2 is the number of quark flavours con¬ 
tributing to the gluon self energy. The collision time 
being of the order of varies from 0.3 — 0.2 fm in 

the temperature range 200 — 400 MeV of QGP. Collision 
time can also be given by -4^ for soft gluon exchange pro¬ 
cesses 0, in which case it lies between 0.5 — 0.25fm for 
the same temperature variation. To ensure the collision 
of heavy quarks with medium particles, r(At) should be 
taken larger than the mean free time rj of heavy quark, 
the latter lying between 1.5 — 0.99fm. The relaxation 
time of heavy quarks is the inverse of the product 
of the gluon de nsity n „(= 16C(3)T 3 /7t 2 ) and scattering 


cross-section a, [17], 19j . Its value is around 1.45 —0.72fm, 


given again the same temperature range. 

The last paragraph contains details of the relevant time 
scales pertaining to heavy quarks. Another crucial time 
scale is the mean free time r^ ed of medium or bath par¬ 
ticles which can be evaluated similarly from the gluon- 
gluon scattering cross section and thereby ranges from 
1 — 0.5fm for the selected temperature variation. To de¬ 
scribe heavy quark transport using FPE, local thermal 
equilibrium of QGP is necessary and to assure it the ob¬ 
servational time scale r should be much larger than rj^ ed . 
Therefore, the following inequality has to be satisfied, 


t^ q > t > rf ed 


( 4 ) 


Only if these time scales are sufficiently well separated, 
the medium effects will be manifested through the trans¬ 
port coefficients. On the contrary, we see that the time 


scales here are of the same order. Thus we are compelled 
to ask whether it is possible to choose an observation 
time scale r within the petite window between and 
. This casts a doubt on the applicability of FPE even 
to a static QGP. Moreover, for an evolving medium, one 
has to ensure temporal local thermal equilibrium. For 
the time being, let us ignore this subtle issue. There is a 
far more serious problem with the accuracy of transport 
equation when applied to evolving QGP as is illustrated 
in the remaining of this Letter. 

Even if the medium is taken to be in temporally lo¬ 
cal thermal equilibrium, the time scale of FPE has to 
be much larger than 0.2fm for QGP. There also remains 
the restriction (mentioned earlier) that the change in 
temperature must not be appreciable during At, whose 
lower limit is set by r. This concept was implicitly used 
in earlier articles fiol . EMI- Let us take r = 0.5fm 
which is midway between and . The corre¬ 

sponding change in temperature in heavy ion collision 
is about 15 — 20MeV [20]. If we calculate the resulting 
uncertainty in time dependent drag and diffusion coef¬ 
ficients (AA, A B) around their respective mean values 
within the temperature range of 200 — 400MeV then as 
per some literature the uncertainty will be not less than 
6 — 10% 13), 21] while others estimate the same to be 


more than 10 — 15% [ll|, |lj] for each. The proportional 


error (A Raa) hr theoretical estimation of Raa due to 
the uncertainty in drag and diffusion values can be eval¬ 
uated from the final momentum distribution ff(jp). For 
an initial momentum distribution centred around po , the 
final distribution is, 
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The above is a solution to the one dimensional FPE 111]. 
The proportional error that Raa acquires is 


A Raa = 


A A 


A 


1 A B 

2 ~B 


5?’ Po (AA,AH) + e 2 (AA, AH, t). 

( 6 ) 

Here, ei contributes to the error through its dependence 
on initial and final momenta of heavy quarks whereas e 2 
contains error accrued due to time dependence. The first 
two terms account for the uncertainty irrespective of the 
momenta of heavy quarks and time. Consequently, we 
find the minimum percentage error or uncertainty in the 
prediction of Raa to be more than 10 — 18% depending 
on different temperature dependences of transport coef¬ 
ficients quoted in recent literature!! [ll-li| HI- Further¬ 
more, this error propagates in time. Thus we see that 


2 We have quoted those articles which together cover the entire 
range of uncertainty. The error calculated from unquoted litera¬ 
ture lies within the specified range. 
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the accumulated error in nuclear modification factor for 
an evolving QGP persisting for 10 — 12fm can be huge 
enough to question the accuracy of FPE. 

To sum things up, we realise that FPE can be applied 
to a slowly evolving medium if its time scale of evolu¬ 
tion is much larger than r. In that case, the temperature 
change is negligible within At and hence, the transport 
coefficients remain unchanged. Then, no question of un¬ 
certainty in Raa arises and we can happily employ the 
clever technique of introducing time dependent functions 
A(t), B(t). But in heavy ion collisions, we are faced with 
the reverse scenario as the medium evolves very rapidly. 
The moderate uncertainty gained by the transport coef¬ 
ficients and the ensuing considerable uncertainty in nu¬ 
clear suppression factor Raa bear testimony to this fact. 
It remains to ascertain whether the accuracy to which 
such transport equations predict heavy quark energy loss 
is satisfactory relative to ongoing experiments with far 
improved precision. Probably no! 

Our argument makes it clear that to mitigate the in¬ 
compatibility of various time scales, the simplest option 
is to think of a description where the time scale r can be 
made arbitrarily small (so that change in temperature 
within that interval is negligible). One such approach 
has been recently devised in the context of heavy quark 
bound states (bottomonium) [22| by modelling them as 
open quantum systems. Development along this line for 
open heavy quarks is deferred for our future investiga¬ 
tion. 


We are thankful to Nicolas Borghini for his critical 
reading of the manuscript. We acknowledge Surasree 
Mazumder and Md. Younus for discussions. T. B. ac¬ 
knowledges support from the DST sponsored project No. 
SR/MF/PS-01/2014- IITI(G) entitled ”A Large Ion Col¬ 
lider Experiment (ALICE) upgrade, operation and uti¬ 
lization” for support. 


* E-mail nirupamdu@gmail.com 

[1] N. Itoh Prog. Theor. Phys. 44, 291(1970) 

[2] J. C. Collins and M. J. Perry Phys. Rev. Lett. 34, 1353 
(1975) 

[3] N. Cabibbo and G. Parisi Phys. Lett. B 59, 67 (1975) 

[4] ATLAS collaboration Phys. Rev. Lett., 105, 252303 

( 2010 ) 

[5] Kohsuke Yagi, Tetsuo Hatsuda and Yasuo Miake Quark 
gluon plasma. From big bang to little bang, Cambridge 
Monographs on Particle Physics and Cosmology (Cam¬ 
bridge University Press 2008) 

[6] E. Iancu, QCD in heavy ion collisions. arXiv: 1205.0579, 
hep-ph (2012) 

[7] F. Gelis. Color Glass Condensate and Glasma. arXiv: 
1211.3327, hep-ph (2012) 

[8] J. D. Bjorken, Fermilab Report No. Fermilab-Pub-82/59- 
THY, 1982 and Erratum (unpublished) 

[9] N. G. van Kampen, Stochastic Processes in Physics 
and Chemistry, North-Holland Personal Library (Else¬ 
vier Science, 2011) 

[10] S. Mazumder, T. Bhattacharyya, Jan-e Alam, S. K. Das, 
Phys. Rev. C 84, 044901 (2011) 

[11] B. Svetitsky, Phys. Rev. D 37, 2484 (1988) 

[12] H. V. Hees, R. Rapp, Phys.Rev. C 71, 034907 (2005) 

[13] S. K. Das, J. -e Alam, P. Mohanty Phys. Rev. C 82, 
014908 (2010) 

[14] Min He, V. Hees, Pol B. Gossiaux, Rainer J. Fries, R. 
Rapp, Phys.Rev. E 88, 032138 (2013) 

[15] M. G. Mustafa, Phys. Rev. C 72, 014905 (2005) 

[16] H. V. Hees, M. Mannarelli, V. Greco, R. Rapp Eur. Phys. 
Jour. C 61, 799 (2009) 

[17] M. H. Thoma, Phys. Rev. D 49, 451 (1994) 

[18] A. Beraudo, A. De Pace, W. Alberico, and A. Molinari, 
Nucl. Phys. A 831, 59-90 (2009) 

[19] T. S. Biro, et. al. Phys. Rev. C 48, 1275 (1993) 

[20] C. Shen, U. Heinz, P. Huovinen and H. Song, Phys. Rev. 
C 84, 044903 (2011) 

[211 S. Mazumder, T. Bhattacharyya, Jan-e Alam, Phys. Rev. 
D 89, 014002 (2014) 

[22] N. Borghini, C. Gombeaud, Eur. Phys. J. C 72, 2000 

( 2012 ) 



